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ABSTRACT 

We consider atomic system of A— type 3-level coupled to 2 mode fields, and derive an effective 
Maxwell-Bloch equation designed for two photon emission between two lower levels. We find axially 



i-t symmetric, topologically stable soliton solutions made of condensed fields. Immediate implication 
— ■ of soliton formation to radiative neutrino pair emission is to enhance its rate, larger than the usual 
oc A^^ factor of target number dependence, along with another merit of increasing the signal to the 
background photon emission. 
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Introduction Superfluorescence(SF) (also called superradiance) is a phenomenon of cooper- 
ative photon emission triggered by macroscopic growth of atomic polarization and induced field fl\ - 
[3]. The total emission rate at its maximum scales with the number of targets as oc N"^ compared 
to the spontaneous decay rate oc A^. SF has been observed in a variety of target states, ranging 
from gas, solid crystals [4], to Bose-Einstein condensates [5], [6]. A simplest system of SF may be 
described by the Dicke model P of 2 levels related by El transition. Its initial state is excited atomic 
state prepared for instance by pulsed laser irradiation, with and without further triggering field. SF 
occurs, in the case of 3 levels, to Raman process, as well [6j. 

The present work has been initiated by our efforts of finding an enhancement mechanism of 
radiative neutrino pair emission from an atomic metastable state |1) |2) + 7 + UiUj (z/j neutrino 
mass eigenstates) [7], [8], in similar ways to SF. This process is useful to determine all neutrino 
parameters, 3 masses and 3 mixing angles, and furthermore to distinguish the Majorana neutrino 
from the Dirac neutrino. After the smallest neutrino mass is experimentally determined, one may 
proceed to detection of relic neutrino of 1.9 K using the Pauli blocking effect [9]. 

After derivation of effective Maxwell-Bloch (MB) equation for related two photon emission pro- 
cess, |1) |2) + 7 -|- 7, we realized that non-trivial soliton solutions do exist, their stability being 
assured by the topological winding number associated with SO (2) rotation of fields. With forma- 
tion of solitons radiative neutrino pair emission is enhanced by O(10^)A^^ for targets such as noble 
gas atoms implanted in solid matrix, opening a new way to perform neutrino mass spectroscopy, 
as envisioned in [8]. Moreover, formation of solitons gives an ideal mechanism of suppressing the 
background process of two photon emission when radiative neutrino pair emission is measured. 

We shall report on this finding of soliton, their fundamental and technological applications, which 
may be wide ranging. Examples that immediately come in our mind for applications include the 
memory storage and quantum computing. 

Throughout this work we assume the natural unit, h = 1 and c = 1. 

Effective Maxwell-BIoch equation A standard method of SF analysis is based on the 
semi-classical set of Maxwell-Bloch equations [4J, involving polarization of atomic targets and induced 
electric fields. We consider a three-level atomic system, |1) (initial state), |2) (final state) and |3), 
with energy level relation of A system, es > ei > €2, €2 = taken for convenience. We take into 
account effects of two fields corresponding to transitions, |1) ^ |3) (pump field Sp) and |3) ^ |2) 
(Stokes field Sg)- The interaction Hamiltonian density is Tij = —{diER^i + d2ER^2) /'^ + h-c (we 
assume two El transitions with dipoles (ij, but extention to Ml transitions should be evident). We 
may derive MB set of equations from the equation for the density matrix, p = —i[H,p] with H the 
Hamiltonian, by using commutation relations, [Rij ,Rki] = SjkRu — SuRkj, along with the Maxwell 
equation in medium. We take the continuum limit so that these transition operators Rij are density 
functions of time and space coordinates. 

When numbers are necessary, we consider two types of examples, the case appropriate for two 
photon emission such as neutral Ba atom, and the other case suitable for radiative neutrino pair 
emission such as noble gas atoms. In Ba atom, 3 relevant states are |1) =^D2, |3) =^Pi (6s 6p), 
|2) =^So. In this case di ~ 0.39 x 10~^cm, ^2 ~ 1.9 x 10~^cm, and energy differences are A31 ~ 
0.83eV, A32 ~ 2.2eV. In the case of noble gas atoms |1) =^P2 , |3) =^Pi , |2) =^So using the LS 
coupling scheme (for instance, the more precise configuration for Xe is |1) = 5p^(^P3/2)6s^[3/2]2 , |3) = 
5p^(^P3/2)6s2[3/2]i , |2) = 5/^5*0), and Xe example gives ^2 ~ 2.2 x IQ-^cm, A31 ~ 0.12eV, A32 ~ 
8.4eV. A precise value of di for Xe is not known, but we may infer it of order a typical Ml transition, 
0[e/2me]. 

We first make an ansatz for field components Ei propagating along z-direction, E^ + iEy = 
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i (£*e-*^f(*±^) + £,e*^^(*±^)) + (h.c.) and for polarization R31 = i?3i,pe-*^f(*±^) , R32 = R32 ,se''^'^*^'^ 
Unconventional sign mixture icUj ,i = s ,p in phases here is chosen for convenience of discussing two 
photon process. 

The original Bloch equation for the matter system is 



e-'^'%ie'^^'Rsi,p) = -^e-^*f;i?2i + -js;B - -j-R,,^, , (2) 
e-^^*a,(e^^^*i?32,.) = -ye-^*£,i?i2 + y^.C - |i?32,. , (3) 
dtB = -^{2di£;Ru,p + rf2^:i?32,s) - + C + n), (4) 
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a^c = -5R(di£;i?3i,p + 2rf2^;i?32,s) - , (i? + c + n) , (5) 



and the Maxwell equation in medium, —d'^E + V'^E = \ {did^R^i + d2d^R32) , where B = R^^ — R 
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and C = -R33 — R22 are population difference, n = n{x) the local number density of target atoms per 
unit volume, and Ai = €3 — ei + Up and A2 = €3 — ujs- For Raman-like processes Aj are taken as 
detuning parameters, and small. Ki oc dfA^^ are El or Ml decay rates corresponding to |3) — >• 

Description of two photon process |1) |2) + 7 + 7 requires another choice for Aj different 
from the Raman process; Ai = A2 = A ~ Aq = (£3 — ei + 63)72, or ~ Up. Rapidly oscillating 
terms oc e^*^^* are averaged out for a long time behavior of variables, and one may drop these terms. 
We then make slowly varying envelope approximation (SVEA) by dropping terms dtRsi against 
Ai?3i, which amounts to balancing equations expressing R-a^i = 1 , 2 in terms of other quantities; 
i?3i,p = (zA + Ki/2)-^diB£*/2 , i?32,s = («A + K2/2)-^d2C£s/'2, . MB equation thus derived involves 
an effective direct interaction of two photon emission, |1) |2) via pump and Stokes field emission; 
frequency dependence indeed gives a correct combination oc £p£*{B — C), with strength did2. 

After a transient time of order the lifetime of the upper level |3), populations of levels approach 
stationary values, namely time independent solution of equations for B , C, eq.(jl]) and eq.([S]), giving 
B = ~nAl^\£,\'^ / D , C = -nA^al^plV^^ D = Ai^al^pP + ^3i\^s\^ , where di\£i\ -C A is assumed. 
The result implies B + C = —n, namely i?33 = (or ^ n). 

Resulting equations are a closed set for two field amplitudes and R21, 
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The field magnitudes are limited by eqs.([I]) - (JSj), and not by eqs.(l6])- (171). This argument suggests 
the maximal magnitudes of fields; \£s\ < 0[\dtRij/{Rijd2)\] and \£p\ < 0[\dtRij / {Rijdi)\], which is 
later related to the soliton mass Ai by \dtRij/Rij\ ~ A^. 

With A 7^ Aq, these MB equations are also useful for description of radiative neutrino pair 
emission, |1) — > |2) + 7 + i/ji/j, with a photon energy set at €3 — A, when the weak term TCw is added 
to the Hamiltonian density and treated as a small perturbation. 
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Axial symmetry and photonic soliton In an axially symmetric case of laser irradiation 
along z-axis, we may introduce cylindrical coordinates, {z,p,9). Spacetime dependence of fields, 
polarization, population difference Sj oc e*"*-''^ , j = s,p,Rs2 oc e*™"'^ , i?3i oc e*"*''^ , is further 
assumed. Consistency of angular dependence in equations and ([7]) requires rus = —mp = m. The 
requirement of one-valued functions demands that m is an integer. This introduces the topological 
winding number m. The transverse operator is then V2 = dp{pdp) / p— m^/p^ . We call this topological 
object the photonic soliton, in short PS [10]. 

Fields have polarizations, and we may use this fact to classify chiralities of solitons. There are 
two types of non-trivial topology of field polarization; (1) TE mode; this is the case explicitly written 
above, E + iEy such that for instance, the Stokes field Eg is ^e«-'^(*±2)+*»"^ times a function 

dependent on the transverse distance p. (2) TM mode; this is the case in which role of the magnetic 
and the electric field is interchanged from TE mode. B = x E has the similar structure to 
TE, hence Eg = —Ey + iE^. Field polarizations are classified by a set of two opposite numbers 
{nis = m,mp = —m), the first entry for chirality of the Stokes field and the second for chirality of 
the pump. 

We further set up an ansatz to work out solutions of equations, ([6]) and ([7j); a functional form 
F{p,z,t){Gs{p),Gl{p)) for )/e^"^^ with separation term ^^^{p) = -{d^F - dlF)/F taken 

more slowly varying with p than time variation. Thus, N"^{p) = Ai^ — n'^ip) consists of two parts, 
where time variation Ai = idtSg/Sg = —idtSp/Sp (soliton mass), and variation along z-direction 
ft(p), taken to reflect effect of index of refraction K^(p) = u^f + A^^(z/^ — l)n{p)/no, with uq a 
central density. This assumption amounts to a physical picture of taking field condensates coherently 
collaborating to propagate with the same index of refraction. 

Using dimensionless quantities, ^ = Aip, Xm = v^Gs/Agg, Y^„i = V^Gp/A^, and 2- 
component notation ip'^ = {Xm , Y^m), one has 

' + ^^^-^ - - - ^m) = no , ^ ( ""f (9) 



where a diagonal 2x2 matrix Q having Qn = uj'^/Ai'^ , Q22 = ^^/Al'^ , and r] = A31/A32 , = 
(i^no/(4A) , Up = c/2?T.o/(4A) , /(^) = n{C,/A4)/nQ. The density profile /(^) = n{p)/no depends on 
how a dense collection of targets (and fields) is excited. We assume that this profile function has a 
characteristic length scale po, which is essentially the soliton size to be determined dynamically. 

Quantum analogy and soliton profile We shall study the case which has the potential of 
observing radiative neutrino pair emission; noble gas atoms (see below on their large rates) implanted 
with a fraction 10~^ in solid para-H2 matrix (po/10 > the lattice constant of matrix ~ 3.8 x 10~®cm). 
The range of parameters are ap/rj = O[0.8(Ne) — 8(Xe)] x 10~^ and the other asf]"^ much smaller, in 
the Xe example asff /oip = O[10^^^]. 

Analogy to quantum scattering problem is useful here. We first note that effect of the right hand 
side (RHS) of eq.([n]) is described by an effective, non-linear interaction. 
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giving (non-linear) propagation of Stokes and pump fields, along with their mixing. This gives a 
Stokes-pump field mixing with effective strength depending on the field ratio r = \£p/£s\ itself; its 
effective strength varies from nd\/ {2A) for small r to —nd\/(2A) for large r. Since ^2 ^ di usually, 
growth of the field ratio is accelerated by non-linear effect, once it is over a threshold value. We 
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shall mainly consider the case of small field ratio (the case of Stokes field dominance), relevant to 
radiative neutrino pair emission. In this case the linearized approximation is excellent, and two fields 
are essentially decoupled. 

We work out in detail the linearized approximation for [Xi , Yli), using the density profile f{C,) = 
^e~^^^°/^o, suitable for the use of the Bessel laser beam of order 1 as a trigger, since it gives an 
interesting scheme of creating fundamental solitons of chirality ±1. The potential V{C,) = 3/(4^^) — 
(z/^ — then has repulsion at the origin, and for a large parameter (i/^ — attraction at 

intermediate region (and a weak repulsion at infinity). Since for a large soliton mass Ai, the energy 
flu is negligibly small, our problem is essentially reduced to finding out (nearly) zero energy solutions 
which exist for discrete set of ^o, thus forming an eigenvalue problem for this parameter. 

For a crude estimate of the eignevalue of zero energy solution, one may use the WKB formula 
for energy levels Es{C,o), namely 2 J^^ d^^jEs — V{S, , ^o) = 27r^s , {s is an integer) with , z = 1 , 2 
turning points, and set the zero energy condition -Es(^o) = to derive eigenvalues of ^o- We thus 
find eigenvalues approximately given by = ^s, with S,s = s^/vr/ \/2(z/2 — 1) . The number of nodes 
for the zero energy solution is of order, O.Sa/ u"^ — I^q. We illustrate in Figure 1 a numerical solution 
of localized field (Xi , V-i)/v^ ~ {£s/ ,£p/^3i)- We confirmed that the WKB energy formula is 
good for large s > O[10]. 




25 50 75 100 125 150 175 200 



Figure 1: Soliton profile. Stokes and pump fields (black) (assumed real) and the potential (dashed 
red) both in arbitrary units vs = Aip are shown for the refractive index u = \/2, the field ratio 
at the origin r(0) = —0.25 (essentially fixed at all ^), and = 12.43 ■ ■ ■ . The barrier region of the 
potential is expanded by 10^ in the inset. 

Relation between the soliton mass Ai and the soliton size po is roughly Aipo = 0[$,s], with one 
of the eigenvalues. The soliton size can be anywhere between atomic distance (~ Inm in solids) and 
target size, maximally the transverse size of laser irradiated region, but most likely sizes of order the 
wavelength of triggering laser are the main component. For numerical estimate below we assume for 
definiteness A4 = ^s7r/(2po), often taking = 1 for crude estimates. 

One might characterize the photonic soliton by saying that it is a stable concentration of fields, 
possibly much below the field wavelength scale, supported by surrounding dense, excited target 
atoms. 

Implication to radiative neutrino pair emission We discuss neutrino pair emission 
|1) |2) + 7 + UiUj, the pair emission caused by the Hamiltonian density, Hw = QwilRzi + (h.c.) 
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with jl the neutrino pair emission current. Features of a single photon energy spectrum, angular 
distribution and distinction from two photon process have been discussed using a different approach 
in [8j. Here we investigate this problem using MB equations, and examine effect of soliton formation. 

The transition amplitude for radiative neutrino pair emission after soliton formation is governed 
by the following equation written in the interaction picture, 

d,R,, = ^(^njt£.l^e-^(^-^--^^)(*±^) , (11) 
4A 1 + -qr^ 

where Ey^{E^) is the energy going into neutrinos (photons). A large soliton mass of ^ A appears 
in the exponential factor. 

The large time limit for the transition rate per unit volume d\R2i{t)\'^ / dt can be worked out by 
using the wave function of soliton g^i^p^^^i^^s Q^^{-M-Ey,,-E^)z ^ resulting in a similar formula to the Fermi 
golden rule, with a difference of large energy factor M. — Eyy — E^. How the total rate scales with the 
total number of target atoms is as follows. We may take the maximal magnitude Eg of order }A/d2-, 
along with r = 0. The enhancement factor per unit atom is K = \£s/ e\^{y /2t['^){2A/ M) ~ 
l/(7rrf2poA3) , where e ~ v^A/W ~ is the field magnitude of a single atom transition, and 

the inverse of the last factor M./2/S. is the number of events for the soliton decay. One has the entire 
enhancement factor for target atoms given by KN"^, where 

,nm, ,eV, n , 10"'^cm, ^ , , 

K ~ 2.4 X 2 _ ^2 

Po A d2 

What happens is formation of soliton of mass Ai, which macroscopically decays. This massive 
soliton goes into many radiative pairs. The microscopic description of this phenomenon is that the 
rate of elementary radiative neutrino pair emission is enhanced by an extra factor, 0[Ai/{d2A^) ~ 
V('^2PoA^)] in addition to the usual coherence factor A^^. 

The elementary rate 7 of radiative neutrino pair emission is roughly given by 
7^,, = G'2,Af2(732/A32)(Ai2/A3i)V(157r5) ~ 3.3 x 10-^''s-\A,2/eVf{^s2/As2)iA,2/A,,)\ For 
noble gas atoms implanted with a fraction 10"'^ in solid para-H2 matrix, the maximum rate is 
4 X 10^^'^s^^A^^ for Ar and 2 x lO^^^s^^A^^ for Xe, giving the total rate of order (40 - 2)Hz for 
A^ = 10^^ (we took = 1 for this estimate). Other noble gas atoms in solid matrix give similar 
rates, somewhat larger for Ne by O[400] than Xe. Alkaline earth and other atoms often give much 
smaller rates. The enhanced rate scales with Af2A|2/(A3iPo) of target parameters, which works to 
give large rates for the A— type noble gas atoms, with large Aj2 and small A31. 

The precise angular distribution of photon in radiative neutrino pair emission depends on how the 
triggering laser irradiation leads to formation of solitons, their number and their size distribution, 
which is a difficult problem to solve. But, the angular distribution from decay of a single soliton can 
be worked out from ffTTl) . Without much calculation we may deduce basic features of photon angular 
distribution by noting combined spatial variation of field and neutrino pair Ssjl- The phase factor 
in the exponent, which needs to be canceled, is ((iA/" + Kp)p + (±k + Kz)z + (m + myy)0)^ with Ki 
the momenta of many neutrino pairs. The correlation to the cylinder axis is evident, and the photon 
emission is confined to a small angle region of ^ < 0[n/M. = (z/^ — l)n(^o)/''^o]- 

On the other hand, two photon process has in RHS of eq. lfTTi) did2n£*£*{\ — 'qr'^)/{AA{l + 
Tfr'^) )e*-^-TT(*='=^) , in which the soliton mass AA is missing in the exponent. Thus, two photon emission 
from solitons do not occur. The two photon process however can occur from amplified pump and 
Stokes field not related to soliton formation, to give a rate simply proportional to A^^. Thus, the 
rate for radiative neutrino pair emission, is more enhanced, at least by the factor r^"^ M. / {d\A^) ~ 
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r~^/(o?2Po^^)5 than for two photon emission, which is of order 10^r~^ or more for noble gas atoms 
in sohd matrix. Incidentally, the elementary rate for two photon emission from a single atom is 
estimated 0[l]sec~^ for noble gas atoms. 

Controled two photon emission, however, becomes possible by using a systematic destruction of 
coherence, such as abrupt modulation of dielectric constant. It should also be noted that stability, 
and possibility of controled coherence breaking, of PS gives an ideal mechanism of enhancing the 
signal to the 2 photon background ratio in measurement of forbidden processes. 

The potential background of multi-photon (more than 2 photon) emission is not enhanced at all 
by soliton formation, thus when the elementary rate of multi-photon QED process is smaller than 
the enhanced rate of radiative neutrino pair emission, say 0[l]Hz, the multi-photon emission does 
not become the major background. 

Applications and outlook Here we briefly discuss some possible technological applications 
using two photon emission caused by controled destruction of stable PS's. 

Topological solitons, both stable or unstable (the case of resonance), are likely to be created in 
the region of high dielectric constant e. One may use for preparation the photonic crystal type of 
medium [12] doped by target atoms of long-lived A— type level such as Ba D-levels. Many photonic 
solitons of small size may be created when an array disk made of rectangular shaped high e material 
is irradiated by the Bessel laser beam of chirality 1 for the trigger. This might serve for the memory 
storage. On the other hand, when a cylinder made of many high e tubes is irradiated by a Bessel 
beam of large aperture, one may expect creation of many PS's of long size, which might serve for 
efficient light transportation. From energetic reasons we expect that PS's of size of order the laser 
wavelength are more likely to be created at its formation. 

Correlation of emitted two pulsed lights after controled coherence breaking is excellent, in direc- 
tion, energy, and chirahty. Correlated emission of strong light pulses after PS destruction may thus 
be useful for quantum information. 

What is pressing is experimental confirmation of the basic idea in the present work, and is to 
clarify how easy or how difficult it is to creat many PS's using material technologically available at 
present. On theoretical side calculation of dynamical time evolution is left to further work. 
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